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Abstract. We state and discuss numerous mathematical identities involving Jacobi elliptic functions 
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1. Introduction 



The Jacobi elliptic functions sn{x,m), cn{x,m) and dn(x,m) with real elliptic modulus parameter 
m (0 < m < 1) have been extensively studied and used in mathematics, science and engineering ||l], ^. 
Recently, while studying ||^, |8| the properties of quantum mechanical periodic potentials [Q, we 
have discovered numerous mathematical identities involving Jacobi elliptic functions. The purpose of 
this paper is to tabulate, derive and discuss these identities. 

To the best of our knowledge, our results are not discussed in the mathematics literature. However, 
we did find that geometrical constructions called the "poristic polygons of Poncelet" give rise to a few 
of our very simplest identities like eqs. (^) and (^) involving just the Jacobi elliptic functions dn(x, m). 
For a discussion of this geometrical approach, see references ||5|, ^. 

Our new identities play a crucial role in obtaining a large class of novel periodic solutions of the 
Korteweg-de Vries (KdV) and modified Korteweg-de Vries equations Q, the nonlinear Schrodinger 
and KP equations, the sine-Gordon equation, as well as the Xcp^ model Q. The solutions obtained for 
the KdV equation all correspond to one gap periodic potentials. This process can be generalized 
to obtain new solvable periodic potentials with a finite number of band gaps |^]. 

If K{m) denotes the complete elliptic integral of the first kind, the elliptic functions sn (x, m) and 
cn(x,m) have real periods 4i^(m), whereas dn(a;,m) has a period 2K{m). The m = limit gives 
i^(0) = 7r/2 and trigonometric functions: sn(j;, 0) = sinx, cn(j;, 0) = cosx, dn(x, 0) = 1. The m — > 1 
limit gives K{\) oo and hyperbolic functions: sn(j;, 1) — > tanhx, cn(x, 1) sechx, dn(x, 1) 
sechx. Therefore, our new identities for Jacobi elliptic functions can be thought of as generalizations 
to arbitrary m of identities involving trigonometric and hyperbolic functions. 

2. Description of the identities 

In all the identities discussed in this paper, the arguments of the Jacobi functions are separated by 
either 2K{m)/p or 4:K(m)/p, where p is an integer (p > 2) depending on whether the left hand side 
of the identity is a periodic function of period 2K{m)/p or AK{m)/p. For any given choice of p, we 
define the quantities Sj, q and di as follows: 

. 2{i-l)K{m) , . 2{i-l)K{m) . , ^ . 2{i-l)K{m) . , . 

Sj = snx + -i — —,m] ,Ci = cn x + ^ — —,m] , di = dn\x + — — —,m] . (1) 

p p p 
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Similarly, we define 

Si = sn[x + — — ^ — -,m] ; a = cn[x + — — ^ — -,m] , ck = dn[x + — — ^ — -,m] . (2) 

p p p 

Each p-point identity which we discuss will involve a cyclic homogeneous polynomial of degree r 
(in Jacobi elliptic functions with p equally spaced arguments) expressed as a linear combination of 
other cyclic homogeneous polynomials of degree r — 2n, where 1 < n < |. We designate this to be a 
p-point identity of rank r. 

Let us consider a few examples to clarify the terminology and establish the notation. A simple 
4-point identity of rank 2 is 

did2 + c.p. = did2 + ^2(^3 + d-^di + d4,di = A , (3) 

where we have used the notation "+ c.p." to denote cyclic permutations of the indices 1,2,..., p. 
Later, we have also used the notation "— c.p." to denote cyclic permutations with alternating positive 
and negative signs. The quantities in eq. (^) are 

di = dn(x,m) , d2 = dn{x + K{m)/2,m) , d-s = dn{x + K {m) , m) , d^ = dn{x + 3K{m)/2,m) . (4) 

Setting X = 0, the constant A can be computed to be A = 2t(l + t^) where 

t = dn{K{m)/2,m) = {l-my^^. (5) 

Similarly, two examples of 3-point identities of rank 2 and rank 3 are 



c\C2 + C.p. = - ^2 ' cid2a3 + C.p. = -g' (ci + C2 + C3) , (6) 



where 

q = dn(2K(m)/3,m) , (7) 

and the arguments are x, x + 4Er(m)/3 and x + 8K{m) /3 respectively. Many more examples are given 
in Tables 1, 2 and 3. Identities of rank 2 are given in Table 1, identities of rank 3 are given in Table 
2 and some examples of identities of rank 4 or greater are displayed in Table 3. 

Although ^-independent constants like A do depend on the number of points p, the rank r, the 
modulus parameter m, and the specific identity involved, for simplicity, we do not usually exhibit 
these dependences explicitly. In fact, the symbols A, B and C appearing in the identities given in 
Tables 1, 2 and 3 are just meant to denote generic constants. They do not all have the same values. 



For illustrative purposes, we now outline the proof of the p-point identity 

did2 + c.p. = A, (8) 

which is true for any integer value p > 1. The left hand side of this identity contains p terms. The 
proof for p = 2 is trivial, since it is well-known that did2 = dn{x,m) dn(x + K{m),m) = \/l — m 
lH, 0]. For p = 3, one needs to compute did2 + c?2'^3 + ^^3*^1 = dn(x, m) dn{x + 2K{m)/3, m) + dn(x + 
2K{m) /3^ra) dn(x + 4i^(m)/3, m) + dn(x + 4i^'(m)/3,m) dn(x,m). This can be accomplished by 
algebraic simplification after using the addition theorem Q 

dn(n + v) = (dnn dnv — msnu cnu snv cnv)/{l — m sn^u sn^v) . (9) 

The result is the constant A for the p = 3 case. One gets A = q{q + 2), where q has been defined in 
eq. ®. 

Similarly, the result for p = 4 has already been discussed following eq. (|3|). In principle, an analo- 
gous algebraic procedure can be used for any value of p, but the algebra becomes increasingly lengthier. 
We have therefore verified identity (|8|) numerically using the advanced mathematical software package 
Maple. Note that for any chosen value of p, the constant A equals p in the limit m = and vanishes 
for m — > 1. 

Given any p-point identity of rank r, one way of generating a new p-point identity of rank r -|- 1 is 
by differentiation and use of the well-known formulas 

— sn(x,m) = cn(x,m) dn(x,m) , — — cn(x,m) = — sn(x,m) dn(x,m) , 
dx dx 

— — dn(x,m) = — msn(x,m) cn(x,m). (10) 
dx 

For example, differentiation of the p-point rank 2 identity (|8|) yields the rank 3 identity 

sici{d2 + dp) + c.p. = Q , (11) 

which reduces to the well-known trigonometric identity 

, 2(i-l)7r, , , 

^sin[2x + -^ ^]=0 (12) 

i=i P 

in the limit m = 0. 

Another p-point rank r identity of interest for r < p is 

did2 . . . dr + c.p. = A (r even) , (13) 



did2 ...dr + c.p. = B^di (r odd) . (14) 

i 

One also has similar identities involving Sj or Cj (instead of di) for any odd value of p. All these identities 
have the remarkable property of reducing the degree of the polynomial in the Jacobi functions from r 
to (1) depending on whether r is even (odd). For small values of p and r, the constants A, B in eqs. 



1^) and (14) are easily evaluated. Some results are 



A{p = 2,r = 2) =2Vl-m , A{p = = 2) = q{q + 2) , A{p = 4:,r = 2) = 2t{l + f) , 



B{p = 3,r = 3) = 3{- 2-1) ' Bip = i,r = 3) = VT 

1 — 



m 



A{p = 4, r = 4) = 4(1 - m) . (15) 
For the special limiting cases m = and m = 1, one gets 

A{m = 0,p,r) = p , B{m = 0,p, r) = 1 , A{m = l,p, r) = B{m = l,p, r) = . (16) 

Another way of obtaining additional identities is by manipulating established identities. For ex- 
ample, for p = 3 and r = 3, eq. (|l4|) is 

did2d3 = j{di + d2 + d3) , (17) 



where B = B{p = 3, r = 3) is as given by eq. ([TsD. Squaring identity (^) for p = 3 and using eqs. (|T5|) 
and ([l^) yields the new identity 

djdl + c.p. = -2(-^ - 1) E dl + [(1 - q'f + - 3 - 4m)] . (18) 

1 — ^ 1 — 

A similar identity is also true for any p, and in fact we have used it in a crucial manner for obtaining 
new periodic solutions of the KdV equation 1^]. However, to establish this p-point identity, one needs 
a generalization of identities dH) and ([l^) . The generalized identities are 

didn + c.p. = A, (n = 2,3,4,...) , (19) 

and 

p 

didj^dj^ + c.p. = Bj2di, (1< ji < 32 < p) , (20) 

i=l 

which we have verified to be true both algebraically and numerically using Maple for many specific 
choices of the integers n,ji,j2,p. 



3. Discussion and comments 



By the techniques described in the previous section, we have obtained a large number of new identities, 
many of which are displayed in Tables 1, 2 and 3. It should be noted that the modulus parameter m 
is not transformed and remains unchanged in all identities. Although it is not easy to give a complete 
systematic classification, we can comment on some general properties. 

(i) For any identity of rank r, the left hand side is a cyclic homogeneous polynomial expression of 
degree r with p terms. 

(ii) If the polynomial on the left hand side is periodic with period 2K(m)/p [4:K(m)/p], then the 
identity involves arguments spaced by 2K{m)/p [4K{m)/p]. 

(iii) The right hand side involves polynomials of rank r — 2, r — 4,... which are "irreducible", some ex- 
amples being J2 di, J2 J2 ^ii J2 ^idii J2 ^i^i, J2 ^^i^ii Yl, CiSidi, etc. and all these irreducibles multiplied 
by df"' where n = 1,2, .... 

(iv) In general, many of the identities of higher rank can be obtained from those of lower rank by either 
differentiation or algebraic manipulation. Similarly, many of the identities of a given rank r (r > 2) 
can be derived from lower rank identities as well as a few identities of the same rank. For example, 
for p = 3, using the identities of rank 2 and three of the rank 3 identities as given by eqs. (|3^ ) to ( |3^ 
one can obtain all other identities of rank 3 as given in Table 2. 

(v) The generic constants A, B, C in any identity can be determined by choosing specific, convenient 
values of x in the arguments. The value x = is a good choice in many cases. Note that for p < 4, 
we have given explicit values for all the constants appearing in the identities - for 3-point identities, 
all constants are expressed in terms of g = dn(2A'(m)/3, m) , and for 4-point identities, all constants 
are expressed in terms of t = dn(i^(m)/2, m) = (1 — m)^/'^ . In writing the constants, we have made 
frequent use of the relationship + 2q^ + (m — l){2q + 1) = 0. 

(vi) Some identities for even values of p involve alternating positive and negative signs. The symbol 
"— c.p." in these identities refers to cyclic permutations with alternating signs. Many of these identities, 
like d1{d2 + dp) — d2{d3 + di) ■ ■ • — dp{di + dp-i) = A{di—d2 + - • - — dp), play a crucial role in determining 
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band edge wave functions of solvable quantum mechanical periodic potentials §. For p = 4 this 
identity is easily derived by starting from identity (|3|) with A = 2t{l + t"^) and t = (1 — m)^^'^. On 
multiplying both sides of this identity by (di — d2 + — 1^4) and using the relations did-s = ^2^4 = 
\/l — m, we immediately obtain the 4-point identity 

diidi + ^4) - c.p. = 2t{l + t + t'^){di - c.p.) . (21) 

(vii) It should be noted that our identities involve cyclic permutations itc.p. of terms which have no 
clockwise or anticlockwise "handedness". For example, for even p, there is no identity of the type 
{did2 — c.p.) proportional to {di — c.p.), since the term did2 has a clockwise handedness. It is only 
when one adds on an anticlockwise handed term dfdp that the the combination [df (^2 + dp) — c.p.] is 
indeed proportional to {di — c.p.). 

(viii) In the limit m ^ 0, one recovers many known non-trivial trigonometric identities. In the limit 
m — > 1, since the period K(l) 00, one usually gets trivial hyperbolic function identities. Both these 
limits serve as a useful check on all the new identities involving Jacobi elliptic functions obtained in 
this paper. Of course, as mentioned previously, software packages like Maple or Mathematica quickly 
provide confirmation of any identity to typically eight digit accuracy. 

(ix) Identities for a given value of p, contain identities of the factors of p as special cases. For example, 
for even p, only half of di, ...,dp are independent since they satisfy identities 

didp+2 = ... = dpdp = Vl — rn , (22) 
2 2 

coming from p = 2. Similarly, the full list of p = 6 identities contains p = 2,3 identities. For example, 
did^ = ^2^5 = dsdQ = \J\ — m and similarly ^1^3 + ^3(^5 + ^5^1 = ^2^4 + di^d^ + ^6(^2 = + 2g, where 
q is as given by eq. (0). 

(x) It should be noted that in many applications like finding new solutions of the KdV equation , the 
identities needed involve summations over all combinations of many (say two) indices i, j = \,. . . ,p. 
These combinations correspond to the sum of several cyclic identities discussed in the tables. 

(xi) In this paper, we have concentrated our attention on cyclic identities in which the arguments 
are separated by fractions of the periods 2K{m) or 4i^(m) on the real axis. However, each one of 
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our identities can be easily translated into a corresponding one in which the arguments are separated 
by fractions of the periods i2K' or i4K' on the imaginary axis, where K' = K{1 — m). For exam- 
ple, the simple 2-point identity dn(x,m)dn(3; + K{m),m) = \/l — m translates to the new identity 
sn(n, m) sn(u + iK' , m) = 1/ ^/m. The general procedure consists of first replacing m by 1 — m [which 



in alternative standard notation Q corresponds to replacing k = ^/m by k' = \J\ — m and K{m) by 
K'], then using the well-known results 

, , , i cniix + K,m) . , dn(ix + K, m) 

sn[x, I — m) = —j^=^= , cn(x, 1 — m 



\J\ — m sn(ix + K, m) ' ' a/I — m sn{ix + K, m) 

dn(x, 1 — m) 



sn{ix + K, m) ' 
and finally changing to a new variable u = ix + K. 

In conclusion, even though Jacobi elliptic functions have been studied for approximately two cen- 
turies, it is exciting to discover new cyclic identities connecting them. What makes our results doubly 
exciting is that the identities play a vital role in the study of periodic potentials § and in yielding 
new solutions of nonlinear differential equations of physical interest Q . 
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Table 1: Identities of rank 2. The symbols A in eqs. ( ^71) and (29) are used generically to denote 
constants independent of x; the constants are in general all different. 

p = 2: 

did2 = Vl-m (23) 



p = 3: [q = dn{2K{m)/3,m)] 



did2 + c.p. = q{q + 2) , ciC2 + c.p. = —— — , S1S2 + c.p. = —{q -1) (24) 

(1 + q)"^ m 

ci{d2 + ds) + c-P- = si{d2 + ds) + c.p. = ci(s2 + S3) + c.p. = (25) 

p = 4: [t = dn{K{m)/2,m) = {l-m)^/'^] 



dids = d2di = yjl-m , did2 + c.p. = 2t(l + t^) (26) 

p = Even Integer: 

did2 + c.p. = A , dids + c.p. = A , ... , didt + c.p. = A (27) 
didz^i = d2dE^2 = • • • = dtdp = ^\ — m (28) 

p = Odd Integer: 

did2 + c.p. = A , C1C2 + c.p. = A , S1S2 + c.p. = A 

didp+i + c.p. = A , ciCp+i + c.p. = A , siSp+i + c.p. = A (29) 

2 22 

ci ((^2 + dp) + c.p. = , si ((i2 + dp) + c.p. = , ci (s2 + Sp) + c.p. = 

ci{dp+i + (ip+s) + c.p. = , si{d,j±i + dp+i) + c.p. = , ci{sp+i + Sp+3) + c.p. = (30) 

22 22 22 
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Table 2: Identities of rank 3. The symbols A in eqs. ( p8[) through (pq ) are used generically to 



denote constants independent of x; the constants are in general all different, 
p = 2: 

did2 ± d^di = \/l — m{di ± ^2) , Cisid2 + C2S2di = (31) 
p = 3: [q = dn{2K{m)/3,m)] 

did2d3 = (g' + - 1) (^^ + ^2 + c^s) (32) 
1 — 

C1C2C3 = ^ _ (ci + C2 + C3) (33) 
-1 

S1S2S3 = nih + h + h) (34) 

1 — 

Ci(s2(i3 + 53(^2) + c.p. = (35) 
cid2d3 + c.p. = -g^(ci + C2 + C3) (36) 



si(i2C?3 + c.p. = 5 (S1 + S2 + S3) (38) 

S1C2C3 + c.p. = ^^^-^^(si + S2 + S3) (39) 



mciS2S3 + c.p. = -(1 + g)^(ci + C2 + C3) (37) 
(2g3_^3g2_2g + 3^_3)^^ 

r 

1 

_ ^ ^ 

dlC2C3 + c.p. = —— -^(di +d2+ ds) (40) 

(1 + qy 

md,~S2h + c.p. = tljl^+A±llL^i^d, + d2 + d,) (41) 

di{d2C2 + d3C3) + c.p. = 2q{q + l)(ci + C2 + C3) (42) 

msi(s2C2 + S3C3) + c.p. = 2q{q + l)(ci + C2 + C3) (43) 
7.T~ 2(g2 + 2g-m+l) , . ,,,, 

di{d2S2 + d3S3) + c.p. = -^^ ^ Si + S2 + S3 44 



-2g(9 + 2) 
(l+g)(l-<?2) 



Cl(c2S2 + C3S3) + C.p. = ^ ^^^^ 2t(^1 + S2 + S3) (45) 



-2q - - - 

Ci{c2d2 + csds) + c.p. = 7— — -y{di + d2 + d^) (46) 

(1 + qy 

msi{s2d2 + hdi,) + c.p. = . — ^ + m) _^ + ^3) (47) 

1 + q 
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dl{d2 + ds) + c.p. = ^\^^ \ di + d2 + ds) (48) 

ci{c2 + C3) + c.p. = ^-^—^ -{C1+C2 + C3) (49) 

-2.~^~N^ 2{q^ + q'^ -q + mq + 2m-l) 

msi(s2 + S3) +c.p. = l-q^ (si + 'S2 + S3) (50) 

cisi{d2 + d^) + c.p. = , cidi{s2 + h) + c.p. = , disi{c2 + C3) + c.p. = (51) 

p = 4: [t = dn{K{m) /2, m) = (1 - m)^/''] 

did2ds lb d2d3d4 + ^3^46?! it ^4(^1^2 = VT^-ni (itdi + ^2 ± ^3 + ^4) (52) 

rf?(rf2 + ^4) ± rfi(4 + c^i) + 4(^4 + ci2) ± dlidi + ^3) = 2t(l T i + i^)(cii ±^2 + ^3^ d^) (53) 

d^ds ± ^2^4 + dldi ± d\d2 = Vl-m {di ± ^2 + ^3 ± d^) (54) 

cisi{d2 + d^) + c.p. = , cisids + c^s^di = , 0282^4 + £454^2 = (55) 
p = Even Integer: 



cisi{d2 + dp) + c.p. = , cisi{d3 + cZp-i) + c.p. = , • • • , cisi(d| + ^2+2) + c-P- = 

ClSl(iE_|_^ + CE^iSE^idl = C2S2<^|+2 + CE_|_2'S£_|_2d2 = • • • = CRStdp + CpSpdt = 

di(d2 + dp) ± c.p. = A((ii ± c.p.) , dl{d3 + dp_i) ± c.p. = A(di ± c.p.) 
di(d| + 6^1+2) c.p. = A{di ± c.p.) , didz^i ± c.p. = Vl — m((ii ± c.p.) 
did2dz lb c.p. = ^(di lb c.p.) , didjdk ± c.p. = ^(di ib c.p.) 

p = Odd Integer: 

For indices 1 < ji < j2 < P '■ 

didj^dj2 + c.p. = A{di + c.p.) , ciCj-^Cj,^ + c.p. = A(ci + c.p.) , siSj^Sjj + c.p. = A{si + c.p.) 
ci{s2dp + Spd2) + c.p. = ci{sp+idp+3 + Sp+sd p+i) + c.p. = 

2 2 2 2 

cid2dp + c.p. = A{ci + c.p.) , • • • , cidp+idp+3 + c.p. = yl(ci + c.p.) 

2 2 

di(d2C2 + dpCp) + c.p. = ^(ci + c.p.) , • • • , di(d£+iCp+i + dp+sCp+s) + c.p. = ^(ci + c.p.) 

2 2 2 2 
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(56) 

(57) 
(58) 
(59) 
(60) 



(61) 

(62) 
(63) 
(64) 



dl(d2 + dp) + c.p. = A(di + c.p.) , • • • , dl(dp+i + dp+a) + c.p. = A(di + c.p.) (65) 

2 2 

cisi((i2 + dp) + c.p. = , • • • , cisi(dp+i + dp+s) + c.p. = (66) 

2 2 

Note that additional identities can be obtained by changing the pair {c,d) in eq. ( |63| ) or eq. (64) 
into any of the pairs (c, s), (s, d), {s, c), {d, c), (d, s). Likewise, additional identities can be obtained by 
changing d to c or s in eq. (|65|) and by changing (c, s, d) to (c, J, s) or {d, s, c) in eq. (|6^) . 
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Table 3: Some identities of rank 4 and above. The symbols A,B,C in eqs. ( |7^ ) through (79) 
are used generically to denote constants independent of x; the constants are in general all different. 



r = 4, p = 2: 

dld2 ± dldi = Vl-m (dj ± d^) , djdl = l-m (67) 
mciSiC2S2 = \/l - m[l - - s^] , cidiC2d2 = -(1 - m)siS2 , sidiS2d2 = -ciC2 (68) 
r = 4, p = 3: [q = dn(2i^(m)/3, m)] 

SiCid2d3 + C.p. = ~^ — ^2 ^("5lCl + c.p.) (69) 

dl{d2 + (is) + c.p. = ^^{dl + c.p.) - 2(1 - m) (70) 
1 — (7^ 

Sl'ilC2C3 + C.p. = ^(Sidi + C.p.) (71) 

1 — 

cidiS2S3 + c.p. = - — ^(ciJi + c.p.) (72) 
1 — (7^ 

"i^ciSiC2S2 + C.p. = ^^^{dj + c.p.) + [m - (2 - m)(l + q)'^] (73) 

r = 5, p = 3: 

o — Imq 

di{s2C2 + S3C3) + c.p. = o (sicidi + c.p.) (74) 

1 — 

~4.~ 2(g2+rn-l) 3 ^ , 2(g3 + g2 ^ mg - g + 2m - 1) 

msi(s2 + S3) + c.p. = {si + c.p.) H 2T2 ('^1 + ^-P-l (^^) 

1 — g (1 — g^)^ 

r = 6, p = 6: 

dl{dld^ + d^ds) + c.p. = A(4 + c.p.) + B{dl + c.p.) + C (76) 

r = 8, p = 6: 

Ci(iiC2d2S3S4S5S6 + C.p. = A{siS2S^SiS5SQ) (77) 

r, p: 

mPs\sl ■■■si = A{sl + c.p.) + B {r = 2p , p = even) (78) 
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For indices 1 < ii < ^2 < • • • < jr-i < P '■ 

didj^dj^ ■ ■ ■ dj^_^ + c.p. = A (r = even , p) 

didj^dj^ ■ ■ ■ dj^^^ + c.p. = B{di + c.p.) (r = odd , p) 
^i^h^j2 ■ ■ ■ ^jr-i + c-P- = A , hcj^Cj^ ■ ■ ■ Cj^_^ + c.p. = A (r = even , p = odd) 
siSj^Sj^ ■ ■ ■ Sj^._^^ + c.p. = B{si + c.p.) , ciCj^Cj^ ■ ■ ■ Cj^_^ + c.p. = B{ci + c.p.) (r = odd , p = odd) 
d\-^{d2 + dp) + c.p. = A{d\-^ + c.p.) + B{d\-^ + c.p.) + ... . (79) 
Eq. ( [79| ) ends with a constant if r is odd and with a term proportional to [di + c.p.), if r is even. 
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